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We show how one can solve the problem of discriminating between qubit states. We use the 
quantum state discrimination duality theorem and the Bloch sphere representation of qubits which 
allows for an easy geometric and analytical representation of the optimal guessing strategies. 
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A fundamental task in quantum information theory is 
that of quantum state discrimination (QSD). It is the 
problem of distinguishing between the possible states of 
a physical system. Generally, state discrimination is not 
a task that can be achieved perfectly; even classically, 
the bias of a die cannot be determined with certainty 
by throwing it once. Moreover, a characteristic feature 
of quantum mechanics is the impossibility to distinguish 
non-orthogonal pure states with certainty. Thus, one 
problem of quantum state discrimination is to determine 
which measurement strategies maximize the probability 
-Pguess of correctly guessing the state. Such measurement 
strategies will be called optimal strategies. Other fig- 
ures of merit have been used to measure the success in 
identifying a quantum state, for example one can con- 
sider the task of unambiguously distinguishing among 
quantum states. However, it is the maximum guessing 
probability P gU ess, and the so-called min-entropy equal to 
— logPguess, which has appeared as the essential quantity 
which measures uncertainty about classical information 
in quantum cryptography, see e. g. 

Since the problem of quantum state discrimination is 
so fundamental to quantum physics, many results have 
been obtained over the last 40 years. The most important 
of these is the formulation of the optimal guessing strate- 
gies given a set ofprobable quantum states, as a semi- 
definite program [3 [| (see also 0]). The semi-definite 
program formulation shows that the quantum state dis- 
crimination problem is computationally tractable. The 
formulation also dashes any hope for an analytical solu- 
tion to the quantum state discrimination problem in gen- 
eral. The special cases which have allowed for particular 
or closed-form answers so far include the well-known case 
of two states (H[ and the scenario where the states them- 
selves can be used to compose the POVM elements of an 
optimal measurement @. Another important subset of 
problems constitute those families of quantum states for 
whom Belavkin's square-root measurement (which was 
later reinvented as the pretty good measurement or 
Belavkin's weighted s quar e-root measurement @,@] is op- 
timal. We refer to [3l. Tirj| - [l3j | and Q for various results 
on the use of the (weighted) square-root measurement. 

The last set of solution strategies pertains to the partic- 



ular problem of distinguishing between qubit states. Ref. 
[lH analyzes the task of optimally distinguishing between 
pure (linearly-dependent) qubit states and finds partial 
results using the Bloch sphere representation. Other par- 
tial results on distinguishing pure or equiprobable qubit 
states were obtained in 16 - We refer to for a cur- 
rent overview of the quantum state discrimination prob- 
lem. 

What is lacking in the qubit state discrimination prob- 
lem is a systematic investigation of how the necessary and 
sufficient conditions of the quantum state discrimination 
duality theorem can be employed to determine the op- 
timal strategies for general sets of qubit states. In this 
Letter we provide such simple analysis. For qubits the 
semi-definite program of QSD reduces to a program for 
finding the minimal enclosing ball of a set of balls in R 3 . 
For this well-studied problem, there exists an algorithm 
whose running time is linear in the number of states to 
be distinguished. Such linear-time algorithm should be 
contrasted with the polynomial but not necessarily very 
efficient complexity of a semi-definite program. In addi- 
tion, it is not hard to show that for qubits there is an 
optimal strategy which involves distinguishing at most a 
subset of four states. We will describe an entirely analyt- 
ical procedure for determining the optimal strategy for 
discriminating between four states. 

We start by recalling the quantum state discrimina- 
tion duality theorem. Let a state discrimination prob- 
lem (X,H,1Z) be defined by a finite alphabet X with 
\X\ = m, an n-dimcnsional Hilbert space % and a set 
of positive semi-definite operators on this Hilbert space, 
K = {Px>0) xex such that Tr [Y, xe x Px] = 1 - For 
such a problem, a guessing strategy is defined by a 
measurement £ = (E x > 0) xeX with Y^xeX E x = In- 
The probability of correctly guessing x using the POVM 
£ is given by p(£) = Y. x ex Tr [E x Px]- Alternatively, 
we can associate the operator a(£) — J2 X E x p x with 
a strategy £ such that the guessing probability equals 
p(£) = Tr [u{£ )]. Hence, the problem of determining the 
optimal guessing strategies is to find an operator a and 
its associated measurement strategies which has maximal 
trace Pgucss • 

Theorem 1 (QSD Duality Theorem [311). Let 
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(X,H,1Z = (px) X £x) be a, discrimination problem. Then 
all optimal measurement strategies £ define the same op- 
erator <j{£) and this operator is called the Lagrange op- 
erator of the problem and denoted as ofo. It has the 
following properties: 



Vx e X : an> p x , 



Vx e X, E x aiz = E x p x <^> £ = (E x ) xeX is optimal 



(1) 

(2) 



(3) 



Alternatively, the Lagrange operator a-n can be found as 
the solution of the following semi- definite program 



min Tr [a] , 
such that Vx € X, a > p x , 



(4) 



and a = . 



Lastly, any guessing strategy £ with operator o~(£) 
Y] E x p x such that 



V.t e X, ° > p x , 



(5) 



is an optimal strategy. 



We will first consider the semi-definite program in the 
qubit setting. Every Hcrmitian operator on qubits can 
be written as p — g {plu + r ■ a) where (p, r) G R 4 with 
r = Tr [pa] and p = Tr [p\. Here 3 = (a x ,a yi a z ) are the 
Pauli matrices. A special feature of 2 x 2 Hermitian ma- 
trices is that the condition of positive semi-dcfinitcncss 
can be expressed as a quadratic constraint, i.e. p > and 
|r| 2 < p 2 . Hence, in this representation the Lagrange op- 
erator for distinguishing the set (p x ) xeX = ((p x ,r x )) xeX 
is given by the solution a = (t, s) € R 4 of the following 
program: 



mint, 

such that Vx G X, t > p x , 
and V.t G X, (t - p x ) 2 > \s- r x \ 2 . 



(6) 



This is an example of a conic quadratic program for vec- 
tors in R 4 (conic quadratic programs form a subclass 
of general semi-definite programs and include linear pro- 
grams [22|]). The constraints in such program are speci- 
fied by quadratic cones. With each vector (p x , r x ) we can 
associate a "light" cone C x which consists of all vectors 
(t, s) such that (t — p x ) 2 > \s — r x \ 2 . This cone con- 
sists of a future light cone C x for which t > p x and past 
light cone C~ for which t < p x . The Lagrange opera- 
tor of a set of states {{p x , r x )) X £x precisely corresponds 
to the point (t, s) which lies in the intersection of the 
future light cones C\ X C X for which the time-coordinate t 
is minimal. Equivalently, the Lagrange operator is the 
state a with minimal trace such that all quantum states 



to be discriminated lie in the past light cone C~ of this 
operator. 

We wish to express the four-dimensional vectors as 
balls in R 3 whose radius depends on the "timc"- 
coordinate which represents the trace of the density ma- 
trix. The mapping will depend on the discrimination 
problem, i. e. let p = ma,x x&x p x > 0. We represent each 
state (p x ,r x ), by a ball B x centered at r x with radius 
\p— p x \. For example, a set of equiprobable states 1Z will 
be mapped onto a set of points in R 3 . Note that each 
ball B x is the intersection of the (future) light cone C x and 
a three-dimensional hypcrplane H p with time-coordinate 
equal to p. Similarly, the ball corresponding to the La- 
grange operator a is obtained by intersecting the past 
light cone of a = (t, s) with the hypcrplane H p . We can 
then characterize the Lagrange operator as follows. 

Proposition 1. Given is a set of balls (B x ) x&x in R 3 
corresponding to a qubit state discrimination problem 
Px = (Px,r x ) and let p = max X £x Px- The Lagrange 
operator a = \{tI-H + s ■ a) of this discrimination prob- 
lem corresponds to the ball B a of minimum radius t — p 
and center S at s which has non-empty intersection with 
each B x , i.e. \fx £ X, B a H B x 7^ 0. This ball with mini- 
mum radius and non-empty intersection will be called the 
"interior ball" of the set {B x ) x ^x- 

The proposition follows readily: since the past light 
cone of the Lagrange operator a should contain all points 
p x , it should have a non-empty intersection with each ball 
B x . The radius of the ball B a obtained by intersecting 
the cone C~ with Tr [a] = t with the hyperplane H p is 
t — p. The Lagrange operator and thus the interior ball 
of (B x ) xe x obtained in this way is always unique. Note 
that for a general set of three-dimensional balls (B X ) X(£X , 
the interior ball is unique if and only if the intersection 
n x B x contains at most one point. 

Before we discuss how to explicitly obtain the Lagrange 
operator and the optimal strategies in the representation 
of Proposition[TJ let us present an observation concerning 
the number of states \X\ = m to be distinguished in a 
Hilbert space of dimension n. If m is large, say m n, 
then it may seem hard to determine the optimal strat- 
egy. In addition, an optimal strategy with m outcomes 
may be hard to implement practically. Fortunately, it 
is not hard to show that there always exists an optimal 
strategy which has at most n 2 outcomes (see also [Hf). 
Davies 2l| has shown that any POVM measurement £ 
on states in a n-dimensional Hilbert space can be written 
as a convex combination of POVM measurements each of 
which has at most n 2 outcomes (this can be viewed as 
a consequence of Caratheodory's theorem applied to the 
n 2 — 1-dimensional real vector space of trace- 1 Hermitian 
matrices). 

Since the guessing probability is a convex (linear) func- 
tion of the guessing strategy, it follows that an optimal 
guessing strategy can be written as a convex combination 
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of optimal guessing strategics each of which has at most 
n 2 outcomes. 

Clearly, the Lagrange operator for this optimal strat- 
egy a = J2xex' E x p x where X' C X and \X'\ < n 2 
obeys Eq. (JSJ) for the ensemble (p x )xex'- Therefore this 
strategy is optimal for distinguishing states in the sub- 
set TZ' = (px)xGX'- Thus aw = &n and all the optimal 
measurement strategies for TZ' are optimal for TZ. 

It follows that in order to find the Lagrange operator 
for the whole set of states TZ, we may consider the La- 
grange operators for subsets TZ' of size at most n 2 . All 
Lagrange operators o-ji' for the subproblems TZ' have the 
property that Tr [(j-r.'} < Tr [on] and thus the Lagrange 
operator o-ji' with the largest trace is the Lagrange oper- 
ator for the whole problem. 

For qubits, these arguments show that one can ob- 
tain the Lagrange operator and the optimal measurement 
strategies by considering the problem of distinguishing 
between at most four states. Hence, if the number of 
states m in TZ is larger than four, one may solve the state 
distinguishability problem by considering all (™) subsets 
1Z\ of four states and compute the Lagrange operator for 
each subset. Among these subsets, one chooses the sub- 
set whose Lagrange operator has maximal trace. This 
Lagrange operator is the Lagrange operator of the whole 
problem and the optimal strategy which discriminates 
the states in this subset is optimal for the whole prob- 
lem. Below, we will present an analytical procedure for 
solving the qubit discrimination problem for four states. 

Note that the procedure of considering all possible sub- 
sets of four elements is not very efficient as a function of 
m. However, one can show that the particular conic pro- 
gram in Eq. ([6]) can be solved by an algorithm which runs 
in a time which is linear in the number of states to be dis- 
tinguished. We note that Eq. ([6|) precisely specifies the 
(unsealed) Lagrange ball centered at s £ R 3 with radius 
t as the ball of minimum radius which includes all (un- 
sealed) balls B x centered at f x with radius p x |25|. For 
this particular problem, i. e. to find the ball of minimum 
radius which includes each of a set of balls in fixed di- 
mension, Ref. p3j presents an algorithm whose running 
time is linear in the number of balls. For the simpler 
case of m equiprobable states represented by points, the 
problem of determining the Lagrange ball is the problem 
of finding the minimal enclosing ball, i. e. the ball with 
minimum radius which contains all points. Finding this 
minimal enclosing ball in R 3 , also known as the 1-center 
problem in R 3 , is an old (26j and well-known problem in 
operations research. It was shown to be solvable with a 
linear time 0(m) algorithm in Ref. [24 1. 

For a general quantum state discrimination problem we 
can distinguish several types of states depending on their 
relation with the Lagrange operator o-ji- We consider 
Eq. ([3]) in Theorem [4] which says that for all a; £ X, 
the POVM elements of an optimal measurement strategy 
are such that E x (ou — p x ) = 0. Since E x > and 



a n — Px > 0, it implies that E x < P x where P x is defined 
as the projector (P x = P x2 ) onto the kernel of an — p x - 

Let us first separate off a trivial case when there exists 
a .To £ A' such that on = p x „ or Kcr(on — p x ) = Iu- In 
this case, not measuring, but simply guessing the state 
to be p Xa is optimal (see also [13] )• Note that this case 
corresponds in the language of qubit cones to the case 
when the forward light cones of p x for x ^ xq all first 
intersect at p Xa ■ If we are not in this trivial case, it is clear 
that for all x, Rank^) < R&nk{P x ) < n - 1. Hence, 
this implies that for qubits, any optimal measurement 
strategy 8 = (E x ) has pure rank 1 POVM elements E x . 
Now we can distinguish several types of states p x : 

There can be states p x for which dim Ker(er-R. — p x ) = 0, 
i. e. E x = P x = 0. We will call such states unguessable, 
because no optimal strategy has X clS el possible outcome. 
For qubits, the geometric interpretation is that for such 
states p x , the Lagrange ball B a is not tangent to the ball 
B x , i.e. the intersection B a D B x contains more than one 
point. An example is the point B in FigQJc). When 
Px is unguessable, we may remove p x from the set TZ 
without changing the Lagrange operator or the optimal 
strategies. 

There can be states for which 1 < dimKer(<7TC — p x ) < 
n — 1. For qubits, this implies that the operator E x is 
proportional to the pure state projector onto this kernel, 
namely E x cx P x = I-u + ^^r^ • a with the Lagrange 
operator a = (t,s). Geometrically, these projectors 
P x are vectors from the center S of the Lagrange ball 
B a towards the center of the ball B x , see e. g. Figs. QJ 
[21 [4] In order to determine an optimal strategy, one 
has to find nonnegative coefficients X x such that for 
E x = X X P X , J2x EX = 7 «- For qubits, the POVM 
completeness condition implies that ~^2 X X x r x = s for 
weights Aa; > 0, J2 X A x = 1. It can be the case that for 
all optimal strategies = 0, hence E x = 0. In this case 
we call p x nearly-guessable; an example is point B in 
Fig. [2(b). It is clear that when p x is nearly-guessable, it 
can be removed from the set of states without changing 
the Lagrange operator or the optimal strategies. In the 
last case when E x ^ for some optimal strategies £ , we 
will call p x guessable. 

Now we present the procedure to discriminate between 
four states. It relies on checking whether the interior ball 
of increasingly larger subsets (from size 1 to 4) intersects 
with the remaining states. If this happens with a subset 
of size k < 4, then there is an optimal measurement 
strategy with these k states as outputs or equivalently 
the k states of the subset are all guessable. W.l.o.g. let 
A be a point and Bb,c,d be three balls. 

1. The point A £ BbC\BcC\Bd in which case the 
interior ball is point A; this is the trivial case where just 
guessing A is one optimal strategy. 

2. The interior ball of a pair of balls, provided that 
it is unique, has a non-empty intersection with the third 
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and the fourth ball. Then, this interior ball is the interior 
ball of the whole set. Clearly, it is sufficient to check for 
pairs of balls such that their interior ball is unique, since 
the Lagrange operator ball itself is unique. Note that 
it is very easy to compute the interior ball of a pair of 
balls by placing the center S of the interior ball midway 
between the two balls. An example of this scenario is 
depicted in Fig. [2] when distinguishing three states. Note 
that in this case, it can happen that the third and fourth 
state are also guessable, for example any two antipodal 
points which are on the blue circle in Fig. [5] would also 
be guessable. If this is not the case, then we consider 
whether: 

3. The interior ball of a triplet of balls, provided that 
it is unique, has a non-empty intersection with the fourth 
ball. Then, again this interior ball of this subset is the 
interior ball of the whole set. We can compute the in- 
terior ball of a triplet of balls, say, Ba, Bb and Be as 
follows. Note that all three balls may have a radius 
strictly larger than zero; if this is so, then we redefine 
the balls as B'a,b,c with radius \p — p x =A,B,c\ where 
p = maji x= A.B,c Px in order to simplify the description 
of the following procedure. It is clear that we can restrict 
ourselves to the plane defined by the centers A, B and 
C and hence we represent the balls by disks T> x= a^b.c 
where w.l.o.g we assume that T> a is a point. The inte- 
rior ball B a can be represented as the interior disk T> a 
(i. e. the disk with minimum radius which intersects all 
disks T> x ). Now, we can first check whether or not the 
interior disk of a subset of two states intersects the third 
state; we already did this computation in Case 2 and 
use the results here. If there is no subset of two states 
whose interior disk includes the third one, then all three 
states are guessable and we determine the interior disk 
as follows: 

We find the center 5 of the Lagrange operator disk 
which is the closest point which is equidistant from the 
disks and the point. Thus, first we find the set of points 
which are equidistant from the point A and the surface 
of the disk T>b , this set of points form a hyperbola Hab 
described by a quadratic equation, see Fig. [3] Similarly, 
we find the points which are equidistant from A and the 
disk T>c, the hyperbola Hac- Computing the intersec- 
tion Hab Pi Hac requires solving a cubic equation and 
the intersection will consists of at most four points. The 
interior disk with smallest radius will have its center S 
at the point closest to A and its radius is \AS\, see the 
example of Figure |U 

Once we have found the Lagrange operator or the La- 
grange ball of a subset of three states, we can check 
whether it includes the fourth state. If no three-state 
subset defines the Lagrange operator, we have as the last 
possibility: 

4. The interior ball of the states is tangent to the 
four balls Ba-, Bb, Be and Bd and we determine the 
center and radius of this ball as follows. The points which 



are equidistant from A and the surface of B x= b (resp. 
B x= c and B x= u) form a hyperboloid of revolution (of two 
sheets) Hab (resp. Hac and Had) in three dimensions 
with the foci A and x = B (resp. x = C and x = 
D). The line-segment of the major axis of Ha x from 
one hyperboloid sheet to the other is of length \r x \ for 
x = B,C, D (see the line-segment of length r in Fig. [3]). 

We can compute the intersection of these three 
(quadratic) hyperboloids Hab, Hac, Had which consists 
of at most 8 points. The computation of the intersection 
is a task which can be reduced to solving two quartic 
equations. Again we choose the point of intersection 
which is closest to A; this is the center S of the Lagrange 
ball and its radius is \AS\. 

Remark: Note that for equiprobable states which can 
be represented as points, the boundary of the minimal 
enclosing ball in Case 2, 3 and 4 is the circumscribed 
sphere of respectively 2, 3 and 4 points. 



A 




(a) (b) (c) 

FIG. 1: S is the center of blue ball B a (here projected on the 
plane defined by the points A, B and C) corresponding to 
the Lagrange operator a. In (a) the triangle ABC is acute- 
angled, and all three states A, B, C are guessable. In (b) the 
triangle is right-angled and B is nearly-guessable. In (c) the 
triangle is obtuse-angled. The interior ball defined by A and 
C includes the point B and B is unguessable. 

In conclusion, we have analyzed the problem of qubit 
state discrimination and have given an analytical proce- 
dure for explicitly solving the problem of distinguishing 
between four arbitrary states. This procedure can be 
easily checked to give the answers derived in Ref. [l6| 
for the problem of distinguishing three mirror-symmetric 
states (see [20|). It is clear that the linear-time algo- 
rithm and the procedure to compute the Lagrange oper- 
ator for a set of four qubit states are particular to qubits 
and do not directly generalize to higher dimensions where 
the constraints of positive semi-definitcness are no longer 
quadratic. However, an interesting question for future re- 
search is to analyze whether the SDP for the Lagrange 
operator, Eq. (U]), is sufficiently simple in its dependence 
on tr and p x to allow for some of the search techniques 
in 0, [24| to apply. In addition, the notion of guessable, 
nearly-guessable and unguessable states and the observa- 
tion about the number of outcomes for an optimal strat- 
egy are general and can also be useful in describing the 
optimal strategies for general quantum state discrimina- 
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FIG. 2: The interior disk of A and T>c intersects the disk 
T>b- The vectors pointing from S to the points A and C' 
represent the pure state vectors E x of the optimal strategy. 
The state ps is unguessable; its disk intersects the interior 
disk of A and T>c and hence E B — 0. Note that the problem 
of distinguishing pA,pB,pc is equivalent to distinguishing the 
equiprobable states represented by the points A,B' and C' . 
Thus the problem of distinguishing non-equiprobable states 
reduces to the problem of distinguishing equiprobable states 
in FigtHc). 




FIG. 3: Let r be the radius of the disk T>b- The point Q is 
the closest point to A which is equidistant to the red circle 
defined by the disk T>b ■ All points equidistant from A and the 
circle form a hyperbola whose foci are A and B, i.e. \AP\ = 
\PT\. In general, the points on a hyperbola are such that 
the difference between the distance to the foci A and B is 
constant. In this case the constant is r. 



tion problems. 

BMT acknowledges support by the DARPA QUEST 
program under contract number HR0011-09-C-0047. 
MED acknowledges support by the department INFRES 
of TELECOM ParisTech. 



* deconinc@telecom-paristech.fr 

* bterhal@gmail.com 

[1] R. Konig, R. Renner, C. Schaffner, IEEE Trans. Inf. Th. , 
vol. 55, no. 9 (2009). R. Renner, PhD Thes is, ETH No. 
16242 (2005), arXiv.org: |quant-ph/0512258| 

[2] A. S. Holevo, Problems of Inf. Trans. 10:4,317-320 
(1974). H. P. Yuen, R. S. Kennedy, M. Lax, IEEE 
Trans. Inf. Theory, Vol. 21, 125-134 (1975). 

[3] Y. C. Eldar, A. Megretski, G. C. Verghese, IEEE 
Trans. Inform. Theory 50, pp. 1198-1207 (2004). 

[4] S. M. Barnett, S. Croke, J. Phys. A 42, 062001 (2009). 




FIG. 4: In this example, the interior disk centered at S is 
tangent to all disks and all states pA, Pb and pc are guess- 
able. The vectors pointing from S to the points A, B' and C" 
represent the pure state vectors E x of the optimal strategy. 



[5] C. W. Helstrom, J. Stat. Phys. 1(2) 231-252 (1969). 
[6] H. P. Yuen, R. S. Kennedy, M. Lax, Proc. IEEE (Lett.), 

vol 58 1170-1173 (1970). 
[7] V. P. Belavkin, Radio Eng. and Electr. Physics, Vol. 20, 

p. 39-47 (1975). 
[8] V. P. Belavkin, Stochastics l,p. 315-345 (1975). 
[9] P. Hausladen, W. K. Wooters, J. Mod. Optics 41, 2385 

(1994). 

[10] A. S. Holevo, "On asymptotically optimal hy- 
pothesis testing in quantum statistics", Theory 
Probab. Appl. 23:2, pp. 411-415 (1979). 

[11] M. Sasaki, K. Kato, M. Izutsu, O. Hirota, Phys. Rev. A 
58 146 (1998). 

[12] M. Ban et al, Int. J. Theor. Phys. 36, 1269(1997). 
Y. C. Eldar, G. D. Forney, Jr., IEEE Trans. Inf. Theory 
47, 858 (2001). S. M. Barnett, Phys. Rev. A 64, 030303 
(2001). 

[13] C. Mochon, Phys. Rev. A 73, 032328 (2006). 

[14] J. Tyson, Phys. Rev. A 79, 032343 (2009). 

[15] B. F. Samsonov, Phys. Rev. A 80, 052305 (2009). 

[16] E. Andersson et al., Phys. Rev. A 65, 052308 (2002). 

[17] K. Hunter, Phys. Rev. A 68,012306 (2003). 

[18] K. Hunter, Results in Optimal Discrimination, Proc. of 
QCMC 2004, Vol. 734, pp. 83-86 (2004). 

[19] S. M. Barnett, S. Croke, Quantum state discrimination, 
Adv. in Optics and Phot. 1(2) 238-278 (2009) 

[20] M. Deconinck, "Optimal Quantum State Discrimina- 
tion", Master's thesis in quantum information theory 
(2010). 

[21] E. B. Davies, Information and quantum measurement, 
IEEE Trans. Inf. Theory, Vol. 24, no. 5, pp. 596599 
(1978). 

[22] A. Ben-Tal, A. S. Nemirovskii, Lectures on Modern Con- 
vex Optimization, Analysis, Algorithms, And Engineer- 
ing Applications, MPS-SIAM Series on Optimization 
(SIAM 2001). 

[23] N. Megiddo, "On the ball spanned by balls", Discrete 
and Computational Geometry 4 (1989) 605610. 

[24] N. Megiddo, "Linear-time algorithms for linear program- 
ming in R 3 and related problems", SIAM Jour, on 
Comp. 12:4, 759776 (1983). 

[25] Note that this representation of the Lagrange operator is 
different (but equivalent in expression) from the interior 
ball representation of Proposition 1. 

[26] In 1857 the mathematician Sylvester posed this problem 
for a set of points in the plane. 



